% 创建时间：2025-06-20 23:15:48
% 作者：兪馳
clc;
clear;

A=[1 2 3;
    4 5 6;
    7 8 0]; % 输入矩阵 A
C=-[10,5,4;
    5,6,7;
    4,7,9]; % 输入矩阵 C
A0=sym(kron(eye(3),A)+kron(A,eye(3))); % 构造系数矩阵
c=C(:); % 将矩阵 C 向量化
x0=-inv(A0)*c; % 求解向量化的解 x0
x=reshape(x0,3,3) % 将向量 x0 恢复为 3×3 矩阵
norm(A*x + x*A' + C) % 验证解的准确性